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' Abstract. This review describes quantum systems of bosonic particles moving on a 

^ . lattice. These models are relevant in statistical physics, and have natural ties with 

probability theory. The general setting is recalled and the main questions about phase 
transitions are addressed. A lattice model with Lennard- Jones potential is studied as 
CO ' an example of a system where first-order phase transitions occur. 

A major interest of bosonic systems is the possibility of displaying a Bose-Einstein 
. condensation. This is discussed in the light of the main existing rigorous result, namely 

. its occurrence in the hard-core boson model. Finally, we consider another approach 

' that involves the lengths of the cycles formed by the particles in the space-time rep- 

. resentation; Bose-Einstein condensation should be related to positive probability of 

CO ' infinite cycles. 

o ■ 



Ph' 1. Introduction 



Statistical Physics is the study of macroscopic properties of systems with a large number 
' of microsopic particles. Its relevance stems from the law of large numbers, allowing the 

state of a system to be specified by the values of a few 'macroscopic variables', although 
^ I the number of microscopic degrees of freedom is enormous. From a probability theory 

point of view, the Ising model of classical spins is an example of identically distributed, 
but not independent, random variables; when couplings are small (high temperature, 
' random variables close to independent), magnetization is zero; for large couplings however 

(strong dependence, or low temperature), the law of large numbers takes a subtler form, 
with two typical values for the magnetization. This behavior is a manifestation of a 
phase transition. Connections between statistical physics and probability theory, such as 
the relation between the physical entropy and the rate function of large deviations, are 



X 



discussed in detail by Pfister in his excellent lectures |^H|. 

While the original motivation for the Ising model resides in quantum mechanics, it is 
considered as a classical model, because energy and observables are functions on the space 
of configurations — in quantum systems, these are operators on the vector space spanned 
by the configurations. There are several reasons for devoting some attention to quantum 
systems. 

• They are closer to the physical reality, and usually of more interest to physicists than 
classical ones. 

• They have richer properties; new types of phases such as superfluidity or supercon- 
ductivity may show up that are intrisically quantum phenomena. 

• They pose a number of mathematically interesting questions. 
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There are three classes of quantum lattice systems. The first class consists of spin 
systems, such as the quantum Heisenberg model, where each site of the lattice hosts a 
spin that interacts with nearest neighbors. In the second class are fermionic systems, 
an example of which is the Hubbard model, where the energy of the quantum particles is 
provided by a discrete Laplacian ('hopping matrix') for the kinetic part, while the potential 
part is given by an operator that is a function of the position operators; particles are 
indistinguishable, so that a permutation of the particles results in the same quantum state, 
up to a sign for odd permutations. The last class consists in bosonic systems that describe 
particles hopping on a lattice and interacting among themselves, but a permutation does 
not alter their wavefunction. There are also other models that have spins and particles, 
particles with spins, or both kinds of particles. 

This review is focussed on bosonic systems. Their great advantage over fermionic ones is 
that they involve only positive numbers, hence natural links with probability theory. They 
also have extremely interesting behavior with various phase transitions, including the Bose- 
Einstein condensation (hereafter denoted BEC), that should be one of the mechanisms 
leading to superfluidity and superconductivity. 

Section ^ introduces the general formalism and defines equilibrium states. This leads to 
the notion of phase transitions, and of symmetry breaking. These ideas are then illustrated 
in a simple boson model with Lennard-Jones potential; its low temperature phase diagram 
is analyzed and shown to display various phase transitions (Section |3|). This can be proven 
by showing the equivalence of this model with a 'contour model' that fits the framework 
of the Pirogov-Sinai theory (Section |^. These techniques, however useful, do not allow 
discussing the occurrence of BEC. We briefly review the main questions in Section ^, and 
state the best result so far — the occurrence of 'off-diagonal long-range order' in the hard- 
core boson lattice model | DLS , KLS |, see Theorem ^T|. We conclude by discussing an 
approach to the BEC that is both geometric and probabilistic, and that involves 'cycles' 
formed by bosonic trajectories in the Feynman-Kac representation. When the temperature 
decreases, the probability of observing an infinite cycle should vary from to a positive 
number, and this transition should be related to BEC. These ideas are described in Section 



2. Mathematical structure 

2.1. Microscopic description. The physical picture is that of a group of bosons on a 
lattice, with the kinetic energy described by a discrete Laplacian, and interacting with a 
two-body potential. 

Let A C Z'^ be a finite volume. The space of 'wave functions' on A is a Hilbert space, 
and a normalized vector describes the state of a quantum particle. For ^' G (8>^=^C^ we 
define the symmetrization operator Sn 

TV 

where the sum is over all permutations of N elements. Then 57v('8'^=iC^) is the Hilbert 
space for bosonic particles, and the Fock space that describes a variable number of 
particles is J"a = e^^QSN{(^^C^)- There is a natural inner product on this space that 
makes it into a Hilbert space. 

This formalism is the natural one from a physical point of view, but it is more practical 
to consider another Hilbert space that is isomorphic to the Fock space above. Thus we 
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start again, this time in the appropriate setting. Standard references are Israel [[sr| and 



Simon |Sim|. 



We consider a Hilbert space Tio; either 7io ~ C°° (more precisely 7io ~ £^(C)), or 
TIq ~ for systems with a 'hard-core condition', i.e. a prescription that sets a maximal 
number N of bosons at a given site. Then we define local Hilbert spaces {'Hx}x£Z<i with 
each TLx — Hq, and for A C Z'^ we set Ha = ^xeA^x- 

A natural basis for Tio is { |no)}noeN; for Ha, an element of this basis is 



\n) = (^xeA\nx), (2.1) 

where n G N^. This represents a state where the site x has Ux bosons. The main operators 
are the creation operator of a boson at site x, noted c^, its adjoint the annihilation operator 
Cx, and the operator number of particles at x, fix = ctcx- Their actions on the above basis 
are 

4 |n) = Vnx + l \n + 6x), 

Cx \n) = ^/n^ \n - 6x) , (2.2) 
fix \n) = Ux \n). 

Here, we denoted \n + 5x) the vector that is equal to ®y,^A \ny + 5xy)- Considering a system 
with hard-core bosons, we demand that c^x\n) = Q \i Ux = N . Notice that the operators fix 
are diagonal in this basis. Without hard-cores, creation and annihilation operators satisfy 
the commutation relations 

[ca;,4] = ^^y (2-3) 

With a hard-core, the relation is 

[cx,cl]=5xy[l-{N + l) (2.4) 

n:nx=N 

In order to avoid extra technicalities associated with unbounded operators, we restrict our 
interest to models with a hard-core condition. 

The energy of the particles is given by an 'interaction', that is, a collection of operators 
H = {Ha)ac1'^ with Ha ■ Ha — *■ Ha- We commit an abuse of notation and still denote Ha 
the operator Ha tX" T^a\a- We define operations {H + H')a = Ha + H'^ and {XH)a = )^Ha, 
and introduce the norm 

= sup^||/7A||e"ll^ll (2.5) 

^ A^x 

for some positive number r, where \\A\\ is the cardinality of the smallest connected set that 
contains A. An interaction is periodic iff there exists a subgroup A' C Z'^ of dimension d, 
such that HrxA = Ha for all x £ A' . Here, Tx is the translation operator. The space of 
periodic interactions with finite norm (|2.5| ) is a Banach space and we denote it B. 
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2.2. Free energy and equilibrium states. The free energ])^ for an interaction H and 
at inverse temperature (3 is 

/(m = -1 lim ^logTr e-^^^cA^^A (2.6) 

where the hmit is taken over a sequence (A„) of volumes such that hm„ ^^jp — '^i 
here, drA = {x G A : dist (x, A'^) ^ r} is an enlarged boundary of A. It is well-known that 
the limit is independent of the way the limit is performed, and that it is a concave 
function of the interactions. 

An equilibrium state pn for the interaction H is a linear, normalized, positive functional 
on the space of interactions, that is tangent to the free energy at H, i.e. for all K £ B, 

PH{K) + f{H)^f{H + K). (2.7) 

To motivate this definition, let us consider the free energy at finite volume f\{H), given 
by (^1^) without taking the limit. The corresponding finite-volume state would be 

oAm-SfF + Ai^l _ ^^'(A^^cAi^A)e-/^^.cA^. 
P^^^^ - dA^^^^ + A=o Tre-Z^S.cA^^A ' 

The definition (|2.7D is therefore more general, and allows to define states directly with 
the free energy in the limit of infinite volumes. The set of tangent functionals at a given 
H is convex; extremal points are the 'pure states'. Existence of more than one tangent 
functionals implies a first-order phase transition. 

A popular definition of equilibrium states in quantum lattice systems involves 'KMS 



states'. They are actually equivalent to tangent functionals, see e.g. [Isi, 3im]. 

One could restrict our interest to operators that are diagonal with respect to the basis 
( |2.1| ) above. In this case, one would consider the configuration space N"^ and the inter- 
actions would be collections of functions on this space. As a result, we have a classical 
system, whose free energy is still given by (p.6|). States can also be defined as tangent 
functionals to the free energy. 

Hamiltonians (or interactions, in our case) may possess symmetries: for instance, a 
translation by a vector of the lattice often does not affect the energy, nor does a rotation 
or a reflection. In quantum statistical physics, one says that U : B ^ B, H ^ H' = U{H) 
is a symmetry if for all volumes A that appear in the limit in (|2.6D there exists a unitary 
operator U\ in TCa such that 

UaY.HaU],'=Y,H'a- (2.8) 

AcA AcA 

Clearly one has f{H') = f{H). 

Let us illustrate this notion on two examples that will be relevant in the sequel. The 
first one is the translation by one site in the direction 1; it is defined by H'j^ = Ha-ci, 
where A — ei = {x : (x(l) -|- l,x(2), . . . ,x{d)) £ A}. Let us assume that the boxes A are 
rectangles with periodic boundary conditions, and 1 ^ x(l) ^ L. Then one can choose 
Ua to be Ua Iua) = |n^), where < = n(^.(i)_i^^(2),...,^.(d)) if x(l) 7^ 1, < = ni^L^^(^2),...,xid)) 
if xfl) = 1. 



^Some authors prefer to define the pressure instead, that is equal to —[3 times the free energy. In 
thermodynamics, the pressure is the potential depending on temperature, volume, and chemical potential. 
It would be physically more appropriate for the discussion of boson models below. The free energy is 
however more convenient for low temperature studies, since lim^^oo f{H) exists in typical situations. 
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The second example is relevant for the Bose-Einstein condensation and is called a 'global 
gauge symmetry'; U\ takes the form U\ = e^"^^GA"^ ^ a G [0, 27r). Hamiltonians describ- 
ing real particles always conserve the total number of particles, and hence possess the 
global gauge symmetry. It can be broken however, yielding states where the number of 
particles fluctuates more than usual.0 We discuss this in Section |5[ 



3. Example: Hopping particles with two-body interactions 

In this section we introduce a simple lattice model and study it by means of geometric 
methods. One obtains that the free energy display angles corresponding to first-order 
phase transitions, see Fig. |2| below. Let us mention that the existence of a first-order 
phase transition in a quantum system in the continuum has been recently established for 
the (quantum) Widom-Rowlinson model |CP, Iof|. 



3.1. The model. The particles have kinetic and potential energy, so that the Hamiltonian 
is 

H = T + V. (3.1) 

The kinetic energy T of particles on a lattice is described by a discrete Laplacian that can 
be written using the creation and annihilation operators in the following way: T = {Ta), 
with 

rp^^\ -t{clcy + clcx) ii A = {x, y} with \x - y\ = I 
[ otherwise. 

We consider here two-body interactions given by a function [/(•) that depends on the 
Euclidean distance between two particles. 



Va 



U {\x - y\) hr^fiy \i A = {x,y} and X ^ y 

iC/(0)ri,(n,-l) iiA = {x} (3.3) 

otherwise. 



The on-site operator \'hx{nx — ^) is the number of pairs of particles at site x, and the energy 
is naturally proportional to it. The model with only on-site interactions was introduced 
in [ FWGF| ] and is usually called the Bose-Hubbard model. 



In order for the Hamiltonian H = T + V to have finite norm ( |2.5| ) , the interaction U 
must have exponential decay for large distances. The density of the system is controlled by 
a term involving a chemical potential, —fiN, where is the 'interaction' that corresponds 
to the number of particles; N^^} = = ii \A\ ^2. 

Let us now discuss in more details the case of a Lennard-Jones type of potential; the 
graph of the corresponding U is depicted in Fig. [l|. We suppose that U{0) = +oo, corre- 
sponding to a hard-core condition that prevents multiple occupancy of the sites. We also 
suppose that the tail 

ur= \U{\y\)\^'^'^ (3-4) 



^ Large deviations of the number of particles in a finite volume are studied in [LLS] in the ideal Bose 
gas, outside the condensation regime. They are indirectly affected by BEC, if the deviated phase is a 
condensate. 
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does not play an important role; only important values of the potential are U{1) and 
C/(\/2)- The results below are valid for Ur ^ uq, the values of uq and r depending on U{1) 
and U{V2). 



U{a) 




a 



Figure 1. The graph of a Lennard- Jones type of potential. 

We start with an analysis of the ground states of the 'classical model' with configuration 
space {0, 1}^ and a Hamiltonian given as a sum over squares S of four nearest-neighbor 
sites: 

HA{n) = J2 [ 2(d -\l ^ n,ny + U{V2) ^ n,.nj [//n, + /i (-l)"nj . 

5CA ^ ' {x,y}cS {x,j/}c5 x&S 



\x-y\=l \x-y\=V2 



(3.5) 



We added a staggered interaction —h{—l)^nx, with (—1)^ = (— l)"^"^. This interaction 
has no physical relevance, but is mathematically useful to uncover the occurrence of phases 
of the chessboard type that breaks the symmetry of translation invariance. One is of course 
interested in what happens when h = 0. 

Four configurations are important, namely (oo)> (oi)' (i o)'^^^ (i i)' respective en- 
ergies are 

g/.,h(10) = f/(^)_/i_ft 



e 



2 2 

'^''^^01) = ;7(V2)-^ + | (3.6) 



e>'^Wl) = 2U{l) + 2U{V2)-f,. 

We make the further assumptions on the potential that U (1) > 0, ensuring a chessboard 
phase to be present, and U{^/2) < 0, so that no phases with quarter density show up 
— they are more difficult to study, since the classical model has an infinite number of 
ground states. In many cases one expects that this degeneracy will be lifted as a result of 
'quantum fluctuations', that is, the effect of a small kinetic energy T. A general theory of 



such effects combined with the Pirogov-Sinai theory can be found in | DFFR , KU |. Notice 
that U{1) > f/(\/2), meaning that at low temperature, the chessboard phase overcomes 
the phase with alternate rows or columns of I's and O's. Energies (|3.6D provide the zero- 



temperature phase diagram and allow guesses for the low temperature situation. 
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3.2. The phase diagram. The situation at high temperature (/? small) is that of bosons 
with weak interactions and no phase transitions may occur. The natural condition for high 
temperature is that /3||i?||r. is small; one can however prove slightly more by not requesting 



that U(0) be small. So we define (compare with (2.5)) 



\H\ 



sup 



.Hl^ll 



(3.7) 



\A\^2 



Theorem 3.1. There exists r < oo such that if P\\H\\* < 1, there is a unique tangent 
functional at H , and the free energy is real analytic in a neighborhood of H . 



This theorem is proven in Section 4.4 using high temperature expansions. We shall 
see below that there may be more than one tangent functionals at low temperature, cor- 
responding to equilibrium states that are not translation invariant. This implies that a 
transition with symmetry breaking takes place when the temperature decreases. Presum- 
ably it is second order, like in the Ising model, but there are no rigorous results to support 
this. 

The limit /3 ^ oo is easily analyzed and is depicted in Fig. |2|. The graph of the function 




Figure 2. The free energy in the limit (3 ^ oo. The phase diagram is 
divided in four domains, corresponding to the empty, chessboards, and full 
configurations. For large (3 and small t, the flat parts bend but the angles 
remain. 

e^'^ is a kind of roof with four flat parts. There are angles between each flat part, so that 
first derivatives have discontinuities there. The two questions that should be asked are: 

• Does this picture survive when adding the tail of the potential, and the kinetic energy 
(hopping matrix)? 

• Does this picture survive at non-zero temperatures? 

The answer to both questions is yes and is provided by the quantum Pirogov-Sinai 
theory. It can be viewed as a considerable extension of the Peierls argument for the Ising 
model. It was proposed by Pirogov and Sinai for classical lattice models |PS, Sin], and 
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extended to quantum models in pKU] , pFE] , |DFFR] , ^Jj\ , Fmjj . These ideas are discussed 
for this model in the next section. One is then led to the phase diagram of Fig. ^. 

Multiple phases and occurrences of first order phase transitions are proven when /3 
is large and t small, i.e. at low temperature and close to the classical limit of vanishing 
hoppings. It is expected that BEG and superfluidity are present in dimension d ^ 3, when 
the temperature is low and with sufficient hoppings |FWGF|. Actually, the situation 
U{0) = oo and U{a) = for a ^ 1 corresponds to the hard-core boson model, when BEG 
is proven at low temperature [ DLS , KLS |; see Section |5|. 



Unicity 



BEC expected 



LRO 



■t 



Figure 3. The phase diagram {t, ^) of the boson model with Lennard- 
Jones potential. There is a unique state (tangent functional) at high tem- 
perature, while a domain with two extremal states, and hence long-range 
order (LRO), is present for low temperature and small hopping (darker 
zone). Most of the phase diagram is not rigorously understood yet. 



The proof of existence of phase transitions were obtained in [ BKU , DFF ] ; it was realized 
in [ FRU that tangent functionals naturally fit in the context of the Pirogov-Sinai theory. 
The zero-temperature energy takes the form (see Fig. ^ 



mm e 

c :) 



(::) 



(3.8) 



where the minimum is taken over the four configurations (qq), (oi)) (lo)' ^^'^ (i i)- 
There are angles at the intersections between different energies. It is not clear whether 
they subsist at finite temperature however — an example where angles disappear is the 
one-dimensional Ising model. The main result of the Pirogov-Sinai theory, in this model, 
is the claim that there exist four functions that are close to the energies ( |3.6| ) , and that 
play the same role: the free energy is given by the minimum of these four functions, and 
hence has angles at their intersections. 

Theorem 3.2 (Free energy at low temperature). Assume d ^ 2. Let U{0) — > oo, U{1) > 
and U{y/2) < 0. There exist Po,r < oo such that if /3 ^ Po and t + Ur ^ 1, there are real 



functions /foV, /f6'l^' /rTi^ ^"^^^ ^^"^ 



■ (U V) 

^0 0' 



■(lO) 



lim 



c :) 



io)'^(i{) 



and 



/3^oo dfi, h (■ ■) dfj,, h 



(::) 
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uniformly in fi,h. Limits are taken in any order. The limit Ur —> means that 
U{a) for all a ^2. 

• The free energy ( |2.6| ) is given by 

c:) ^■■> 

• The functions are in ^, h with uniformly bounded derivatives. Furthermore, f^'^.^ 
is real analytic in fi, h when f^.'^s^ is the unique minimum. 

The phase diagram is therefore governed by these four functions; clearly, it is symmetric 
under the transformation h — > —h. Let be the coexistence point of (qq) and the 
chessboards, i.e. 

J (0 O-s — J (0 l-s^ \'^-^) 

and ^2 be the coexistence between the chessboard and (i {). There are exactly two 
extremal tangent functionals for < < //2 and /i = 0. Exactly three for fj, = fii and 
h = 0, as well as for fi = fi2 and h = 0. There is a unique tangent functional everywhere 
else. 

Among the consequences are various first-order phase transitions. For instance, 



d 



. . (3-10) 



d 



fli,h 



(3.11) 



dh h=o- dh" h=o+ 

for m < n < IJ.2; also, if /i = 0, 
and similarly at ^2- 

Construction of the functions ('metastable free energies' in the Pirogov-Sinai termi- 
nology) is done in two steps. First, using a space-time representation of the model, one 
defines an equivalent contour model. This step is explained in the next section; it gives 
the opportunity to make the link with a stochastic process of classical particles jumping 
on the lattice. The second step is to get an expression for the metastable free energies 
starting from a contour model, and this is achieved using the standard Pirogov-Sinai the- 



ory [PS, Sin]. This is only outlined here. Ideas are described e.g. in | Kot ]; we also mention 



| Uel | for a self-contained review which includes precise statements on tangent functionals. 



3.3. Incompressibility. The space-time contour representation actually allows us to ob- 
tain more. The total number of particles is conserved, and as a consequence the ground 
state of the quantum model has same density as that of the model without hoppings, and 



hence the compressibility is zero. The following observations were made in [BKU2|. 

Since a state is a linear functional on the space of interactions, we have to understand 
what is the density of the systems. We consider the interaction N: 

^ ' (3.12) 
otherwise; 

if p denotes a state, than the corresponding density is p{N). It is a function of the chemical 
potential p. One defines the compressibility kt, 

f^T = ^p{N) (3.13) 
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where the derivative is with constant temperature (i.e. (3). The theorem below claims 
incompressibility of the ground state, and also that the low temperature states are close 
to incompressible. It holds in all dimensions. 

Theorem 3.3. Let U{0) oo, U{1) > and U{y/2) < 0. There exist /?o,r < oo such 
that if (3 f3o and t + Ur ^ 1, one has 

\p{N)-po{N)\^Ce-''''- 

for some C < oo, r' > 0. 



4. The space-time representation and the equivalent contour model 



4.1. Equivalence with a stochastic system. We start with the finite- volume expres- 
sion for the free energy, 



rfl,h 

J A 



-i-logTr e"^^^cA-f^A 
/3|A| 



(4.1) 



with H = T + V — fj,N. Notice that the last two interactions are diagonal with respect to 
the basis (|2.1[) . 

One can give various probabilistic interpretations for ( [4.1[) , see e.g. [ Tothf . A natural 
one is a continous-time Markov chain where the collection of random variables {n{t)}t ^ o 
take values in {0, . . . , N}^. Let us introduce the set of 'neighbors' of a configuration n: 



J\f{n) = {n' : 3a;, y with |x — y| = 1 and n'^ 
The generator of this random process is 



^zx + ^zy for ah z e A}. (4.2) 



1 if n' G AA(n) 

-\Af{n)\ if n' = n 
otherwise. 



(4.3) 



The partition function Z\ = Tr e /^^aca^a ^]^g expectation 

r/3 



= IE[o,/3]('x[«(0) = n(/?)]exp|- / dr^^ U{\x - y\)n^{T)ny{T) - p^n^ir) | 
^ x,yeA xeA 



(4.4) 



Another representation that is more appealing for the physical intuition involves continu- 
ous-time simple random walks. It was explicite d in |CS|| and used to obtain a bound on 
the free energy of the Heisenberg model [ |CS2| , fTothf ! Let {xj( 
random walks with generator 



iit)}t^o, I ^ j ^ N,he 



L 



xy 



1 a \x — y\ = 1 
—2d if X = y 
otherwise. 



(4.5) 
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Then the partition function takes the form 



N >0 



Nl 



XI 



E EE 

-/3 



!^7r(i) 



1 ^ i ^ TV 



exp{-^ dTj2U{\Xi{T)-Xj{T)\)} 



Xi{0) = Xi,l ^ N). (4.6) 



Here particles have to start and end in A, but they are meanwhile free to move outside. 
One could impose more stringent boundary conditions, by defining a generator L^y that 
does not allow particles to leave or enter A, or by adding an infinite potential outside 
of A. It is however useless, as the free energies corresponding to these various partition 
functions have the same thermodynamic limit. 

Notice the sum over permutations in (^^) ; this suggests to consider probability on sets 
of permutations, for instance the probability that the permutation has infinite cycles. We 
discuss this in Section ^, where ( [4.6| ) is heuristically important. 

Let us mention another example of close ties between quantum systems and probability 
theory: Aizenman and Nachtergaele showed the equivalence of a quantum spin chain with 



a stochastic process, which is itself equivalent to a two-dimensional Potts model [AN]. 
Using results established for the latter, the authors can draw new conclusions on the 
former. 

4.2. Equivalence with a contour model. A way to derive these stochastic represen- 
tations is by using Duhamel formula: if A and B are two matrices, then 

Jo 

= e^ + / dri . . . dr^ e'^^ B e^^^-n)^ ^ _b e^^-"™)^ . (4.7) 

„ ^ 1 iO<ri<...<r„<l 

Here we set A = X]aca(^^ ~ /^-^A — ^Pa), with P denoting the staggered interaction, 
and B = Y^aca'^a- Taking the trace, and introducing 1 = J2n right of 

each operator B, we get the following expression: 



m>0 



E E / dTi...dr^e-^^^^("i)(ni|rAjn2 

l,...,Am ni,-,nm -'0<ri<...<r„</3 

where we introduced H\(n) 



^(/3-r„)HA(ni) 



(4.8) 



^^(-^(n| iJ^ |n). One recognizes (4^) and ( |4.6|) . Indeed, 
the sum over {Ai} is actually over pairs of nearest-neighbors; (ni| Ti^n^^y] 1^2) is zero unless 
712 is a 'neighbor' of ni, i.e. it is the same as ni up to one particle that moved from x 
to y, or from y to x. Finally, { e"^'^^"'^^"^^^'^^"'^^ } is represented in (|4.4D and ( [4.6D by the 
exponential. 

To each choice of m, {^j}, {nj}, {Tj}, corresponds a space-time picture illustrated in 
Fig. ^. We write n(T) the configuration at time r, that is, n(T) = Uj if Tj-i ^ t < Tj. 

The goal is to extract some information on the analytic properties of the free energy, 
that is, the logarithm of the partition function. A technique that was proposed in 1975 
for the study of extensions of the Ising model is the Pirogov-Sinai theory [P^ 



was later extended to quantum systems in [BKU, DFF, DFFR, KU, FRU|. The strategy 
is to map the quantum system onto a 'contour model'. The latter is a model where the 
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Figure 4. Space-time representation of the partition function as expanded 
in ([4.8D. There are four bosonic trajectories in this picture. 

states are not configurations or vectors of a Hilbert space, but sets of mutually disjoint 
contours; the statistical weight e~^^ is replaced by a product of individual weights for 
each contour. 

Let us describe in details the setting of a contour model. 

A contour ^ is a pair {A, a), where ^ C Z'^ is a connected set and is the support of A. 
In order to define a, let us introduce the closed unit cell C{x) C M'^ centered at x; the 
boundary B{A) of A C 7,'^ is 

B{A) = {C{x) n C{y) ■.xeA,y^A}. (4.9) 

The boundary -6(^4) decomposes into connected components; each connected component 
b is given a label «{, S {1, . . . ,p}, and a = (ab)- 

Let K (ZTy finite, with periodic boundary conditions. A set of contours {^i, . . . 
is admissible iff 

• Ai G A for all i, and dist {Ai, Aj) ^ 1 if i 7^ j. 

• Labels aj are matching in the following sense. Let W = A \ Uj^^Aj] then each 
connected component of W must have same label on its boundaries. 

For j £ {1, . . . ,p}, let Wj be the union of all connected components of W with labels j 
on their boundaries. 

The partition function of a contour model has the form 

Za= Y1 11^-4,) He-Z^'^^l^'l, (4.10) 
Mi,...,^fc}i=i i=i 

where the sum is over admissible sets of contours in A. 

The weight w{A) of a contour ^ is a complex function of the temperature and of the 
parameters of the phase diagram (here /u and h) that is real anlaytic in all these parameters. 
Furthermore, we need that 

\w{A)\ < e-^^«l^l e-^l^l (4.11) 

for a large enough constant r (depending on d and p). This typically holds when P is 
large. We also need that partial derivatives of the weights with respect to and h satisfy 
the same bounds. 
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Many classical lattice models have such a representation. The usual way to define a 
contour model is to attribute a set of contours to each configuration. One is given a finite 
set of periodic configurations ('low energy configurations', or 'reference configurations'), 
and one defines 'excited sites' as those sites whose neighborhood does not agree with 
any of the reference configurations. The set of excited sites decompose into connected 
components, that are supports of the contours. Outside the contours the configuration 
agrees with one of the reference configurations, and the labels indicate which one. 

The labels are important because the weight of a contour typically depends on which 
configuration lies outside. If we want this weight to depend on the contour only, we need 
to provide the information contained in the labels. 

We are looking for a similar approach here with the space-time representation. On the 
one hand, we expect the phase diagram to display four phases: a phase with very low 
density, corresponding to (go); two chessboard phases, ( J ? ) and ( ? J ); and a phase with 
density close to 1, (} i)- These are our reference configurations. On the other hand, we 
suppose here that particles have small hoppings, so that jumps are typically rare in Fig. 

In order to get contours that have supports on a lattice, we discretize the continuous 
direction. Let /3 such that (i = Mf3 with M an integer. We consider the lattice A = 
A X {1, ... , M}. A site x = {x, s) £ A is 'in the state ( o q)' if for all y with \y — x\ ^ 1, 
and all {s — < r < s/3, we have ny{T) = 0. We make similar definitions for the other 
three reference configurations. 

Cells that are not in such a state are excited. Connected components of the set of excited 
cells are the supports of the contours, and labels take values in {(oo),(oi))(io)'(i i)} 
and contain information on which configuration touches the support. This is illustrated 
in Fig. |5|. 



I I I I I I \ I \ I \ I \ I \ I \ I \ ^ 

Figure 5. Contours in the space-time representation. The contour on the 
left separates the empty configuration from a chessboard one, while the one 
on the right is due to the motion of a particle. 

Summing first over contour configurations, then integrating over compatible space-time 
configurations, we can rewrite ( [4.8D as 

Za= n^A)ne-'^^^'^'l. (4.12) 

{Ai,...,Ak}i=i i=i 
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The expression for the weight w{A) is comphcated, but the exponential bound ( [4.11| ) is 
not too hard to obtain. It will require f3A to be large, and (3t to be smah. Theorem 
is then a result of the Pirogov-Sinai theory, see for instance |Ue]]. 



4.3. Consequences of the contour representation. A few words need to be added in 



view of Theorem 3.2. The density is 

Tr 

p = lim 



Za 



(4.13) 



This expression for the density agrees with that in terms of derivative of the free energy, 
provided the latter is differentiable. Indeed, let /(/i) be the infinite volume free energy as 
a function of the chemical potential. It is concave, and if it is differentiable at n we have 

p=-^ = - hm (4.14) 
d/i A d/i 

The space-time expansion of ( 4.13| ) was studied in | BKU2 |. Due to the conservation of 
the total number of particles, differences between the density of the quantum model (with 
hoppings) and the classical one (without hoppings) lead to contours that wind around the 
torus A X [0, /3]per- Hence their length is at least /3, and no such contours survive when 
taking the limit /3 — > cxd. As a consequence, the density of the quantum model is locked 
to the classical one. 

This clearly implies that the compressibility vanishes at zero temperature. To obtain 
the low temperature bounds requires some more work, that also goes through an expansion 



involving winding contours | BKU2 | 



4.4. Proof of Theorem 3.1. We conclude this section by proving that there is a unique 
equilibrium state at high temperature, as stated in Theorem |3.1| . It strongly relies on ideas 
discussed above, with many simplifications. We show the equivalence between the quantum 
model and a polymer model — this is a contour model without labels (i.e. p = 1). Once we 
have obtained this equivalence, the results follow from cluster expansions [ [KP| , |Dob , BZ]. 
Using the Duhamel formula (4.7), we get 



Tr e-'^^^cA-f^A ^rjVe-Z^E.eA^^w + ^ (_l)"^ ^ f dn . . . dr^ 

m^l Ai,...,A™-^0<^i<-<^-</3 

Q--riT,a: ff^^Q-{-r2-'ri)T,^ H^.^y _ _ g- (/^-Tm ) E a: -'^{^} . (4.15) 



Let Fx^^ be defined by 



(0) 



Tr e"^^{-> , 



(4.16) 



with the trace taken in the single-site Hilbert space n^. We also set fJ^^ = J2xeA We 
define polymers as connected components of the set U'jL^Aj, and the weight of a polymer 
A to be 



w{A) = e<'' ^{-ir E / dn...dr^ 

m^l Ai,...,A„ JO<Ti<...<T^<P 

\Aj\ ^ 2,UjAj=A 

Xr e~'^^^^s^"'^{^> e~^^2~^^-*^^"'^{^> . . . e~''^~'^™-*^^6'4-'^{^} . (4.17) 
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Then 

k 

{Ai,...,A,}j=l 

This is the partition function of a polymer model. 

We need a bound on the weight of the polymers. Since the dimension of the Hilbert 
space Ha is (A^ + 1)1^1, we can estimate the last line of ( |4.17[ ) by 

m 

Tr • ^ (Af + 1)1^1 W II e-'^^f-} || \{ \\Ha, ||. (4.19) 



Furthermore, 



so we obtain 



-^^{^} II ^ e-^^^"' , (4.20) 



K^)|^ (AT + 1)1^1 e-^l^l 11^^^ 



m^l A'cA 



\A'\^2 

^ g|A|log(7V+l) g-r|A| ^\A\\\H\\; ^ ^^ 21) 

This satisfies the assumptions of the cluster expansions when ||-ff||* ^ 1 and r — log(A^ + 
1) — 1 is large enough (depending on d only). One then obtains an exact expression for 
the infinite- volume free energy: in the translation invariant case {Fx = Fy and w{A + x) = 
w{A)), the mean free energy is given by 

k 

/ = ^o°^-^ E ^^{M,...,A,)f{w{A,), (4.22) 

with the sum over clusters, that is, fc-tuples (^i, . . . , A^), k ^ 1, such that their union 
yj^^^Aj is connected. The combinatoric factor ip^ {Ai, . . . , AjS) has an expression involving 
the graph of k vertices with an edge between i and j whenever Ai U Aj is connected. The 
results on cluster expansions include bounds ensuring the convergence of the sum (fl.22|); 



see e.g. |KP, Dob, BZ| for detailed results and proofs. 

By averaging over a cell whose dimensions are given by the periods of the interac- 
tions, one obtains a similar expression in the case of periodicity rather than translation 
invar iance. 

If ll-ff II* < 1, then \\H + \P\\* < 1 for all perturbation P, and A in a neighborhood of 0, 
and one can perform the above expansions. As a result, we obtain a free energy /(A) that 
is given by a convergent sum of clusters, with weights that are analytic in A. Therefore 
/(A) is real analytic, and there is a unique tangent functional at H. 



5. A DISCUSSION OF THE BOSE-ElNSTEIN CONDENSATION 

5.1. The origins. The story started in 1924 when Bose sent a paper to Einstein, that 
was previously rejected by Philosophical Magazine. Einstein translated it into German 
and recommended its publication in Zeitschrift fiir Physik; he wrote articles shortly after- 
wards in Sitzungsherichte der Preussische Akademie der Wissenschaften (1924-25). The 
'Bose-Einstein statistics' for quantum particles (in particular photons) was uncovered, and 
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a curious phase transition was proposed, where the ground state of the one-particle Hamil- 
tonian is macroscopicaUy occupied. This is the Bose-Einstein condensation for the ideal 
boson gas (that is, without interactions). 

For some time it was not clear whether such a transition was really occurring in the 
nature; but London proposed in 1938 that superfluidity in Helium was a consequence of a 
Bose-Einstein condensation, an idea that is largely accepted nowadays. 

Is there a condensation for interacting systems as well, and what does it mean? These 



questions were addressed by Feynman [Fey]; he proposed the idea that the transition 



corresponds to positive probability for the occurrence of infinite cycles in the space-time 
representation — this will be discussed in greater details in the next section. Feynman's 
conclusion is that weakly interacting systems behave like non-interacting ones, albeit with 
a larger effective mass, and still display condensation. 



Direct experimental evidence of BEC has been observed only recently [ AEMWC ] 



5.2. General ideas. A system of bosons in the continuum is described by the Hamil- 
tonian 

H = -^E^^ + E^(\^^-^^\)- (5-1) 

j=l i<j 

Here, Aj is the Laplace operator Yla=i « li and xj is the position of the j-th. particle. 



The low temperatures should be described by the Bogolubov theory, see e.g. [ Lieb , ZB | for 
an introduction and partial justifications. The Bogolubov theory relies on the assumption 
that most of the particles are in the ground state of the Laplace operator (that is, the 
Hamiltonian for the ideal gas), which is false in presence of interactions. Still, many 
predictions are correct; in particular, it gives a value for the ground state energy per 
particle cq at low density, 

^^^27r^ 3 
m 

where p is the density and a is the scattering length of the potential U. This formula 



has been rigorously established by Lieb and Yngvason iLYj. This and other results are 
reviewed in [ Lieb2| . 



Further developments led to the concept of off-diagonal long-range order due to Penrose 
and Onsager |POf| . Take e.g. the lattice model of Section ^. One considers the following 
order parameter: 

Tr c e ^-^cA -'^'^ 

where H = T + V — fiN, and the traces are in the Hilbert space ®x(^k'Hx- Here, it is 
natural to set periodic boundary conditions for A. The question is: 

Does\\m\x-y\^oo{cl:Cy) differ from 0? 

The equilibrium state at high temperature is unique and clustering, see Theorem 3.1 , 
and hence BEC must be searched at low temperatures. 

5.3. The hard-core boson lattice model. There is one rigorous result concerning the 
existence of condensation in a reasonable model of interacting bosons. This is a lattice 
model where bosons interact with hard-core repulsion, i.e. the Hamiltonian ( p.l| ) with 
U{0) — > oo and U{a) — > if a ^ 1. The theorem below is due to Dyson, Lieb and Simon 
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|DLS], and Kennedy, Lieb and Shastry []KLS |. It is stated for 3 or more dimensions and at 



low temperature, but it also holds for the ground state of the 2-dimensional model [KLS] 



Theorem 5.1. Take 3, H = T + V with U{0) oo, U{a) ^ /or a > 0. Then there 
is (3o < oo such that for (3 > I3q, 

lim {clcy) / 0. 

\x-y\^QO 

This theorem implies the existence of a phase transition in the sense that the state 
(•) is not clustering. It is established using 'reflection positivity', introduced in | FSS|| for 



proving spontaneous magnetization in the classical Heisenberg model; its difficult extension 
to quantum systems was done in pLSj . The claims of |DL|, that are relevant here 



deal with spontaneous magnetization in the spin ^ x-y model. Let us discuss analogies 
between spins and hard-core boson systems. For the latter, we take — and define 
self-adjoint operators {S^\s'i'\s'i'^}^^'j^d, that commute if they are located on different 
sites, and satisfy [s'i'\s^^] = iS'i^^ (and permutations of (1,2,3)) at a same site. (These 
matrices are called Pauli matrices.) The x-y model has interaction —Slps'^^ — S^^ si^^ 
on nearest-neighbor sites x, y, and zero otherwise. 

The correspondence to boson models is done by setting 

ct=5«+i5(2) 

= - iSP (5.4) 

nx = + i 

In the case of hard-cores (with = 1) the commutation relations are [cf, c^] = if x 7^ y, 
and {c^,ci} = 1, where {•,•} denotes the anticommutator. It is easy to check that these 
also follow from the commutation relations of spin operators, and from definitions above. 
The x-y model is equivalent to H' , 

^f^^ i-ki'^Uy + clc^] iiA = {x,y}, \x - y\ = I 
1 otherwise. 

Off-diagonal long-range order is then equivalent to spontaneous magnetization in the 1-2 
plane. 

5.4. BEC & symmetry breaking. The Bose-Einstein condensation is related to a sym- 
metry breaking, namely 'global gauge invariance'. Let us note that the Hamiltonian ( |3.lD 
conserves the total number of particles, i.e. 

[^HA,^h,]=0. (5.6) 

AcA xeA 

Therefore one can define the unitary operator U\ = e'^^^e^"^ ^ which is a symmetry of 
the Hamiltonian. Its action on creation and annihilation operators is 

UacIUx'= e'^cl 

Uac^U^' = e-"" cl (5.7) 

This is easily seen from the action of all these operators on elements of the basis ( p.l[ ). 

To study the properties of the free energies as a function of the interactions, one has to 
proceed similarly as in Section ^. Recall that we added a non translation-invariant (and 
non-physical) interaction HP and looked at a phase diagram where /i is a parameter. This 
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is similar here. First, we need an interaction that does not conserve the total number of 
particles. The simplest choice with self-adjoint operators is Q = (Qa), with 



Qa 



e^" cl + e- 







if A = {x} 
otherwise. 



(5.8) 



Supposedly, there is a unique tangent functional to the free energy at H + hQ for all h ^ 0, 
but there should be an infinite number of extremal states at H, if the temperature is low 
enough; each of these extremal states is indexed by a E [0,27r). Since there is a unique 
equilibrium state at high temperature (Theorem 3T), we face here the breakdown of a 
continuous symmetry. It should occur at low temperature and if the dimension of the 
lattice is greater or equal to 3. 

There is no rigorous result to support this discussion, besides the weaker — but impor- 
tant! 



statement of Theorem 5.1 in the case of the hard-core boson gas. 



6. Infinite cycles: context and conjectures 
6.1. Heuristics. In the last section of this brief review, we discuss an approach to the 



BEC initiated by Feynman 50 years ago [ Fey | , that focusses on the occurrence of infinite 
cycles in the space-time representation. Its appeal to probabilists should be evident — it 
looks at first sight like a percolation phenomenon. However, the one-dimensional nature of 
cycles makes them harder to study than clusters. Still, some progress should be possible. 

The partition function for the Hamiltonian ( |5.1| ) can be expanded via Feynman-Kac; 
setting h^/2m = 1, the partition function is given by 



N >0 



Tre^jv i=l Xi(l3)=x^(^i) 

drY,U{\x^{r)-x,{T)\)]. (6.1) 



Here, integrals are over Brownian paths starting at xi and ending at See | |Gin | for an 

introduction to functional integration. This expression is very similar to ( |4.6| ) for lattice 
systems and is illustrated in Fig. ^ 




Figure 6. Feynman-Kac representation of the partition function for 
bosons in the continuum. The picture shows a situations with five par- 
ticles and two cycles, of lengths 4 and 1. 
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The space-time is periodic in the vertical direction, so it is topologically equivalent to a 
cylinder. Bosons wind around the cylinder, forming cycles (see Fig. |6|). Feynman's idea is 
to consider the length of these cycles, and to look at the probability of occurrence of infinite 
ones. He identifies the onset of a positive probability to a Bose-Einstein condensation. In 



his paper [Fey] he argues that interactions only slow down the diffusion of bosons, without 
forbidding infinite cycles, and he concludes that BEC should also occur in interacting 
systems. 



Cycles were studied in [Siito], where it is proved in particular that, in the case of the 



ideal gas (that is, non-interacting particles), infinite cycles do occur below the transition 
temperature for BEC. The converse statement, namely absence of infinite cycles in absence 
of BEC, is not proven yet, although it is doubtlessly true. 

However, the equivalence between BEC and occurrence of infinite cycles is not obvious. 
Consider e.g. the model discussed in Section ^. Our results imply absence of BEC at 
low temperature and with small T; on the other hand, even though they have restricted 
motions, bosons can interchange with neighbors, and infinite cycles seem likely for low 
enough temperature, if the dimension is greater or equal to 3 — this has something to do 
with probabilities of recurrence of random walks. A lattice model can be viewed as a con- 
tinuum model where the particles have condensed (in the usual sense) and are displaying 
long-range order. The following conjecture is compatible with these considerations: 

Conjecture. 

• Occurrence of BEC implies positive probability of infinite cycles. 

• Positive probability of infinite cycles, and absence of long-range order, imply 
occurrence of BEC. 

In the hope of shedding some light on this discussion, we introduce a simple lattice 
model of cycles, state some (rather obvious) properties and propose some conjectures. 



6.2. A simple lattice cycles model. The expression ( |6.lD for the partition function 
starts by an integration over all initial positions of the particles; let us suppose that they 
are located on the sites of the lattice Z*^ — assuming that density fluctuations do not play 
an important role in the onset of BEC, this assumption is a mild one at low temperature. 
Furthermore, we replace the integral over Brownian paths by an effective weight 



n 



e 

(with A C finite). A natural choice for (^{x,y) is \x — y\'^/(3, with | • | the Euclidian 
distance, and /3 the inverse temperature. Indeed, the Brownian paths for a time interval 
[0,(3] diffuse like y^. Other choices are possible, for instance \x — yp / P with 7 > 2 to 
account for large interactions. One could also simplify the problem and consider 



\i X = y 

1/(3 if |x- 2/1 = 1 (6.2) 
00 otherwise. 



In any case, we restrict the choice of ^ to one that satisfies 

< 00, (6.3) 



E 



ensuring that particles do not jump to infinity in one step. 
Let us describe carefully these cycles models. 
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The lattice is Z'^, and we denote by B the set of bijections Z'^ — > . Given x,y G Z'^, 
let Bxy = {vr G IB : 7r(x) = y}; then we define B' to be the algebra made out of all such 
sets and their complements. 

Next we set B(A) = {vr G B : 7r(a;) = x for all x ^ A} the set of permutations that are 
trivial out of A. Since B' is countable, there exists a sequence of boxes A = (A„)„ ^ o such 
that for all B £ B' the following limit exists: 

lim — y I k G si rr e-«(^'"(^» = P{B). (6.4) 
Asa Z(A) ^ 

^ ^ 7reB{A) xeA 

The normalization Z{A) is 

Z{A) = Yl e-«(^'"(^)) . (6.5) 

7reB(A) xeA 

The probability ( |6.4| ) extends to the smallest cr-algebra generated by B' , that we denote 
B. 

A cyc/e is a sequence c = (xi, . . . , X|c|) of different sites; we identify (x2, . . . , xi) = 
(xi, . . . , X|c|). The set of permutations Be = {tt £ M : vr(xj) = Xj_|_i,l ^ j ^ n} (with 
X|c|+i = 2;i) is an element of B, and the set of cycles is countable. Therefore, the set 

B^=M\U B, (6.6) 

is also in the cr-algebra B. It represents the event 'the origin belongs to an infinite cycle', 
and is the central object of our attention. 



6.3. Few results and important conjectures. There are no infinite cycles at high 
temperature; the condition of the following theorem is easy to check for small /3. 



Theorem 6.1. // 



1=1 

En 

c90 j=l 



then P{Boo) = 0. 



Proof. Let i?>n be the set of permutations where the origin belongs to a cycle of length 
greater than n. One has 

5>i D S>2 D . . . and B^o = S>„. (6.7) 

n 

Then P{B^) = lim„P(S>„). Since 

S>„= y U B^, (6.8) 

\w\=n 
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with w = (0, xi, . . . , Xn-i, x) a self-avoiding walk from to x, and Byj = ^^=iBxj_^,xj-, 
one can write 

P(i3>„) = J] JimPA(i3^) 

x^O «!:0— ►x 
|u!|=n 

x^O ui:0— >x 
|ui|=n 

= lim V n e-«(^--^^) 

|c|>n 
\c\>n 

The first inequality is Fatou's lemma. The last term goes to as n — > oo since the sum 
over all cycles containing the origin converges. □ 

The typical picture at high temperature is that of Fig. |^ (a). Most cycles involve a 
unique site and have length 1. When the temperature decreases, cycles lengths should 
increase, as depicted in Fig. |^ (6). The cycles model resemble that of multiple random 




(a) {h) 

Figure 7. Expected typical configurations of cycles, (a) at high tempera- 
ture and {h) at low temperature. 

walks interacting via exclusions. Assume for a moment that ^ is given by (|6.2| ) with (5 = oo, 
that is, cycles have nearest-neighbor jumps. One can generate a configuration of cycles 
by starting at the origin and doing two self-avoiding random walks in different directions. 
When they eventually met, we close this cycle and start another pair of walks from a free 
site, that have to avoid the first one. One repeats the procedure until all the sites have 
been considered. This actually does not give the same probability distribution on the 
configurations of cycles, but one can expect similar behavior. There is a natural question 
in this process: Is there a chance that after n steps the two legs have not crossed? If the 
non-crossing probability remains finite when n goes to infinity, there are infinite cycles. 
It is actually known that the random walk is recurrent in dimension 2 and transcient 
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in dimension 3 and higher. Considerably extrapolating this argument, one obtains an 
illustration on the fact that BEC occurs only in dimensions greater or equal to 3. This 
also suggests the natural conjecture that infinite cycles do occur in this model at low 
temperature and d ^ 3. 

Acknowlegments. I am grateful to R. Moessner and Y. Velenik for a critical reading of 
the manuscript. 
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